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ABSTRACT 

We show that using the redshift dependence of the deep cluster abundance 
to detect the nature of Dark Energy is a serious challenge. We compare the 
expected differences between fiat ACDM models, with different l^mo, with the 
difference between ACDM and dynamical DE models. In the former case, cluster 
abundances in comoving volume and geometrical factors act in the same direction, 
yielding a significant difference between the expected angular densities. On the 
contrary, when we keep a constant and change the DE nature, abundances 
in comoving volume and geometrical factors act in the opposite direction, so that 
the expected differences in angular densities reduce to small factors. 

Subject headings: cosmological parameters 

1. Introduction 

High redshift supernovae, data on the cosmic microwave background (CMB), as well as 
on the large-scale galactic distribution (Riess et al. 1988, Perlmutter et al 1988, Tegmark et 
al. 2001, De Bernardis et al 2000, Hanany et al 2000, Halverson et al 2001, Spergel et al 2003, 
Percival et al. 2002, Efstathiou et al 2002) indicate that ~ 70 % of the world contents are 
due to a smooth component with negative pressure, such that the ratio w = p/ p <~ —0.8. 
The nature of this compenent, dubbed dark energy (DE), is still open for debate. Candidates 
range from false vacuum, yielding a positive cosmological constant A, to a self-interacting 
scalar field (Ratra & Peebles 1988; Wetterich 1988) to even more exotic physics of extra 
dimensions (e.g., Dvali & Turner 2003). 

While the observed value of A implies a dramatic fine-tuning of vacuum, at the latest 
phase transition, ACDM models provide an excellent fit to data. Alternative viable DE 
models must, first of all, rival their success. Therefore, their phaenomenology must be 
and is hardly distinguishable from ACDM. This calls for tests able to discriminate between 
different DE natures. 
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In recent work, the evolution of the cluster mass function has been shown to have a 
significant dependence on DE nature. This has been shown first on the basis of a Press 
& Schechter (PS) formulation (Mainini, Maccio & Bonometto 2003), then using n-body 
simulations (Klypin et al. 2003; see also Linder & Jenkins 2003), which confirmed PS 
findings. 

There are a number of other cosmological measures which can contribute to discriminate 
between different DE natures. This is important also on the light of the point that we wish 
to make in this paper. In fact, while the number density of clusters n{> M, z), in a comoving 
volume, has a significant z dependence, the measurable signal is much less significant. On 
the light of observational uncertainties, it seems unlikely that the nature of DE can be easily 
discriminated along this pattern. 

In a sense, this contradicts the expectations for a number of deep sky surveys (see, e.g., 
Davis, Gerke & Newman 2004) based on the ancient intuition of Hubble (1926), that the 
number-redshift relation can be used to determine the geormetry of the world. Discrimi- 
nating the DE nature is much harder than fixing the matter density parameter Q,^ in a flat 
ACDM model. 

It is not so because of a reduced impact of DE nature on geometry or on n(> M,z), 
but because their effects on geometry and n{> M, z) tend to erase each other. 



Let us consider a family of objects whose (cumulative) mass function is n{> M,z). 
Indipendently of the actual z dependence, the angular number density of objects belonging 
to such family, with redshift between z and z + Az, in a spatially flat geometry, reads 



2. 



Geometrical effects 




(1) 



with D{z) — dr/dz. For the models we are considering, it is useful to show that 




(2) 



and, therefore. 




(3) 
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In fact, for flat models, on the past light cone, adr = —cdt and, therefore, by dividing the 
two sides by dz — —da/a^, we obtain adr/dz — a^cdt/da, so that: 



D(z) — ^ — — ^ fhere, as usual, H = —\ (4) 
dz H[z) \ a/ 



In turn 



dr c 

-- I here, as usual, 
') ^ 



3 Qm{z) 3 ^m{z) 



„^2 ^mo(l + Zy . 



SO that eqs. (2) and (3) follow. 



These equations lead to handable expressions when DE has a state equation w — p/ p 
with constant w. It is then easy to see that 

D{z) = {c/H,){l + z)-^'^[n^o + (1 - ^mo){l + zf^]-''\ (6) 

while r{z) can be worked out by integrating from to z. We can also easily differentiate 
both with respect of Jl^o and with respect to finding that 

^ l-{l + zf^ 

^hJ (1 + ^)3[Q^„ + (1 - + zf-'Y ^ ' 

and 

dP'' ^ c 2 (1 - ^mo)31n(l + (1 + zf^ 

dw ^Ho^ (1 + Zfi^mo + (1 - ^mo){l + zf-^^ ' 

SO that we expect that the geometrical factors increase when Q^deo — 1 ~ ^mo increases, 
whenever w < 0, and decrease when w increases. Notice that eq. (8) is not a functional 
derivative, as the expression (6) holds just for constant w. 

An extension to dynamical DE can be however performed by using the interpolating ex- 
pressions yielding ^^(2;), for RP (Ratra & Peebles, 1988, 1995) and SUGRA (Brax & Martin 
1999, 2001, Brax, Martin & Riazuelo 2000) models, provided by Mainini et al (2004). In 
these cases DE is due to a scalar field, self-interacting through a potential whose expression 
depends on an energy scale A (see Appendix A for details on the dynamical DE models 
considered here). Using the interpolating expressions or, equivalently, direct numerical inte- 
gration, we obtain the results shown in Figure 1. Here the z dependence of the geometrical 
factor D{z)r'^{z) is shown for three ACDM models (ACDM 06, ACDM 07, ACDM 08, with 
^mo = 0.4, 0.3 and 0.2, respectively), as is obtainable from eq. (6). In the same Figure 
we show the z dependence of geometrical factors also for SUGRA and RP models, with 
^mo — 0.3 and with the A parameter fixed at lO^GeV. Ho is 70km/s/Mpc in all models. 
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In the absence of number density evolution, Fig. 1 would also show the dependence of 
the angular number density on z. 

The ACDM models considered are characterized by quite different values of Q^o- Among 
them, only 0.7 approaches available data, which are coherent with an interval not wider 
than 0.02-0.03. In spite of that, the discrepancy between ACDM 07 and ACDM 06 only 

marginally exceeds the discrepancy between SUGRA and ACDM 07 and is smaller than the 
difference between RP and ACDM 07. Notice, in particular, that the RP model can be 
considered as discrepant from data as ACDM 06 or ACDM 08. 

Altogether, Fig. 1 shows how strongly the nature of DE affects geometrical factors and 
that, in the presence of a non-evolving population, an insight into the DE equation of state 
can be provided by the z dependence of their angular number density. 



3. Cluster number evolution 

A fair insight into the evolution of the number of clusters with the redshift z^ can be 
obtained by using a PS expression. Sheth & Tormen (1999, 2002) as well as Jenkins et 
al (2001) provided expressions more closely fitting n-body simulations (which are already 
reasonably approached by PS results). The latter expressions are more complex and include 
more parameters, while their use is unessential for the present aims. 

The expected (differential) cluster number density n{M) , at a given time, is then given 
by the expression 

M 

j{v)vd\ozv = —n{M)Md\o%M . (9) 

Pm 

Here is the matter density, v — 5c/ a m is the so-called bias factor, M is the mass scale 
considered, um is the r.m.s. density fluctuation on the scale M and 5c is the amplitude that, 
in the linear theory, fluctuations should have in order that, assuming spherical evolution, full 
recoUapse is attained exactly at the time considered (in a standard CDM model this value 
is ~1.68; the difference, in other model, ranges around a few percent). As usual, we took a 
Gaussian f{y) distribution. 

Together with eq. (9), we must take into account the virialization condition, which yields 
significantly different density contrasts Ac in different DE models. Further details can be 
found in Mainini et al (2003). 

In Figure 2 we show the cumulative cluster number density, n{> M,z), obtained by 
integrating n{M), for the same models of Fig. 1 , for M = IO^^Mq/i"^. All models are 
normalized to the same cluster number today and the redshift dependence of n(> M, z) 
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is clearly understandable, on qualitative bases: When ACDM models are considered, the 
evolution is faster as we approach standard CDM. Hence, ACDM 08 (ACDM 06), being the 
most distant (nearest) model to standard CDM, has the slowest (fastest) evolution. ACDM 
07 just stands in between. 

When we compare dynamical DE models to ACDM 07, we know that they also yield a 
slower evolution, i.e. structures form earlier in models with w > —1. 

Accordingly, apart of numerical details, the behaviors shown in Fig. 2 reasonably fit our 
expectations. 

4. Conclusions 

Let us now put together the results shown by Figs. 1 & 2 and evaluate the angular 
number densities of clusters in redshift intervals z, z + Az with Az = 0.1. If ACDM models 
are considered, the geometrical factor D{z) r^{z) and the evolutionary factor n{> M, z) vary 
in the same directions, when is modified. 

On the contrary, for dynamical DE models, the two terms vary in the opposite directions 
and, therefore, model differences tend to erase. 

Accordingly, the results shown in Figure 3 are essentially expected. It is then easy 
to argue that testing the differences between cluster numbers is really challenging, when 
different DE models are compared. On the contrary, different flmo^^ cause comparatively 
huge shifts. Data on cluster masses are obtainable either from temperatures T or luminosities 
L (see, e.g., Pierpaoli et al 2004). Finding a fair cluster mass is perhaps the main reason 
for uncertainty, but also the cluster redshift determination is subject to errors. In Figure 4 
we compare the expected mass functions, for the ACDM 07 and the two dynamical DE 
models considered, at two mass scales {W^MqH-^ and 4.2 ■ W^M^h^^), showing also the 
uncertainty in the halo mass function caused by a 5 % indetermination in the mass: 6N = 
[Mn{M)/N{> M,z,Az)] (SM/M) , still for Az = 0.1. Of course, 6M/M ~ 0.05 is not 
easily reachable, on the basis of the present capacity to reconstruct cluster dynamics. The 
plots also assume that a complete cluster sample in a solid angle of 1 steradian is available. 

Using deep cluster distributions, to discriminate DE models, is therefore hard. A com- 
parison between the two panels of Figure 4 however shows that differences become wider 
as we go to smaller masses. Accordingly, on group mass scales, a further enhancement can 
be expected. An even stronger difference can be expected, when angular number densities 
of galaxies are being compared. In particular, comparing the two panels shows that, when 
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lowering the mass scale considered, differences in trends (vs. redshift) are enhanced. On 

galaxy mass scales, therefore, the very peaks of distributions could fall at differend redshifts. 
Predictions for galaxy scales cannot be formulated only using a PS recipe, but should include 
an accurate consideration of how large halos split into galactic objects. This will be the topic 
of future research. 

Before concluding let us however outline that, here above, the dependence on z of 
cluster features has been assumed not to depend on DE nature. As a matter of fact, when 
wc consider different state equations for DE, we change the z-t relation and the same redshift 
corresponds to a different time. This is true also when is changed, of course. However, 
while, in the latter case, this change of the z-t relation has a counterpart in a change of halo 
angular densities, there is no equivalent counterpart in the former case. Comparing the z-t 
relation, obtained from a study of cluster morphologies or galactic evolution, with number 
density predictions can be therefore a way to discriminate between different DE natures. 
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A. Dynamical DE models 

Dynamical DE is to be ascribed to a scalar field, 0, self-interacting through an effective 
potential whose dynamics is set by the Lagrangian density: 

IIde --\^9 {d^W + V{cl>)) . (Al) 

Here g is the determinant of the metric tensor g^j^^, = a?{T)dXndXi, (r is the conformal time). 
In this work we need to consider just a spatially homogeneous (f) {dicf) <^ <j); i — 1,2, 3; dots 
denote differentiation with respect to r); the equation of motion then reads: 

+ 2^0 + a2^ = O. (A2) 
a acp 

Energy density and pressure, obtained from the energy-momentum tensor T^,^, are: 

, = = £ + ,^lTt = f-^-V(^), (A3) 

SO that the state parameter 

P 'ti^/2a + V('t>) 
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changes with time and is negative as soon as the potential term 1^(0) takes large enough 
values. 

The evolution of dynamical DE depends on details of the effective potential V{4>). Here 
we referred to models proposed by Ratra & Peebles (RP: 1988, 1995), yielding a rather slow 
evolution of w, and Brax & Martin (SUGRA: 1999, 2001, see also Brax, Martin & Riazuelo 
2000) yielding a much faster evolving w. Altogether, RP and SUGRA potentials cover a 
large spectrum of evolving w. They read 

V{<P) = RP, (A5) 

V{(j)) = —— exp(47rG02) SUGRA. (A6) 

These potentials allow tracker solutions, yielding the same \ow-z behavior, almost indepen- 
dently of initial conditions. In eqs. (A5) and (A6) , A is an energy scale, currently set in 

the range 10^-10^° GeV, relevant for the physics of fundamental interactions. The potentials 
depend also on the exponent a. Fixing A and a, the DE density parameter Qde is deter- 
mined. Here we rather use A and ^2^^ as independent parameters. In particular, numerical 
results are given for A = lO'^GeV. 

The RP model with such A value is in slight disagreement with low-/ multipoles of 
the CMB anisotropy spectrum data. Agreement may be recovered with smaller A's, which 
however loose significance in particle physics. The SUGRA model considered here, on the 
contrary, is in fair agreement with all available data. 
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Fig. 1. — Redshift dependence of geometrical factors for ACDM and dynamical DE models. 
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Fig. 2. — Cluster number, in comoving volumes, in ACDM and dynamical DE models. In 
this and in the following plots A2; = 0.1. 
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Fig. 3. — Cluster number, on the celestial sphere, in ACDM and dynamical DE models. 
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Fig. 4. — Cluster number, on the celestial sphere, in ACDM and dynamical DE models with 
^mo = 0.7. Effects due to an uncertainty of 5% in the mass determination are also shown. In 
the range 1 < z <2 and for M = lO^^/i~^M0, discrimination may be possible after suitable 
improvements of data and theory 



